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^ . Abstract 

The density of states and energy spectrum of the gluon radiation are calculated for the 
^ I color current of an expanding hydrodynamic skyrmion in the quark gluon plasma with a 

^> ' semiclassical method. Results are compared with those in literatures. 
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1 Introduction 



In this letter, we address the issue of gluon radiation during the hydrodynamic stage in the 
evolution of the deconfined hot QCD matter or quark gluon plasma (QGP) [1] (for review 
see for example [2]). 

The medium induced gluon radiation has been thoroughly explored in the context of final 
state partonic energy loss or "jet quenching" [3]. The spatially extended nuclear matter 
affects the processes of fragmentation and hadronization of the hard partons produced in the 
relativistic heavy ion collisions. Essentially all high p± hadronic observables are affected at 
collider energies and the degree of the medium modification can give a characterization of the 
hot QCD matter in the deconfined phase. In principle, the medium induced radiation effect 
emerges from thermal QCD per se. However, in practice, different approximation schemes 
are applied giving consistent results [4, 5]. On the other hand, gluon radiation has also been 
considered in the context of gluon density saturation in the initial stage, where a strongly 
interacting gluonic atmosphere is crucial for the rapid local thermalization for the deconfined 
QCD matter [6]. 

The time evolution of the RHIC "fireball" can influence the observable particle production 
spectra. Given a strong initial interaction, the resulting state of matter is usually modeled as 
a relativistic fluid undergoing a hydrodynamic flow. Generalized fluid mechanics that charac- 
terizes the long-distance physics of the transport of color charges has been developed for this 
purpose [7] (for review see [8]). Recently, we discovered a type of single skyrmion solutions in 
color fluid [9]. Moreover, we found an interesting case in which the time-dependent skyrmion 
expands in time, which is in accordance with the expanding nature of the fireball generated 
in RHIC experiments [10]. The pattern of gluon radiation pertaining to the color current of 
these non-static configurations is an important character of this color skyrmion. So in this 
letter we calculate this radiation spectrum in a semiclassical approach. The main results from 
our calculation are the following. There is a fast fall-off in the UV side of the spectrum but a 
smooth peak dominates the intermediate energy. And in IR, a long tail is the characteristic 
feature. 

The organization of this paper is the following. In Sect. 2, after a brief review of the 
nonabelian fluid mechanics, we calculate the nonabelian current corresponding to the soliton 
solution. In Sect. 3, semiclassical gluon radiation is calculated. In Sect. 4, comparison of the 
radiation spectrum in our hydrodynamic approach and in other approaches is carried out. 
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2 Color current of an expanding soliton 



Given the thermalization of hot QCD matter above the deconfmement transition temperature, 
the transport of the color charges in the volume of the nuclear size can be modeled by a 
nonlinear sigma model in a first-order formalism 

C = fu tl -F{n)- 9eff J a »Al. (1) 

This nonlinear sigma model describes an ideal fluid system. The configuration of this fluid 
is described by a group element field U, which shows up in the velocity field 00 M 

w„ = ~Tr{<r 3 U%U). (2) 

Conjugate to the velocity is the abelian charge current . It is easy to see that the first term 
in the lagrangian density (1) gives rise to the canonical structure of the fluid system. The 
fact that we will consider only one abelian charge current means that U takes value in an 
SU(2) group. The information about the equation of state (EOS) of the fluid is contained in 
the second term, which is essentially the free energy density of the fluid. In fact, energy and 
pressure densities are given by the ideal fluid formula 

e = F, p = nF'-F. (3) 

Here n is the invariant length of , n 2 = j^j^ . The third term is the gauge coupling of the 
fluid with an external gluon field A a ^ with an effective coupling g e ff- J afl is the nonabelian 
charge current which is related to the abelian current by the Eckart factorization J afl = Q a j' J ' 
where Q a is the nonabelian charge density of the fluid configuration 

Q a = ^Tr(a 3 U^a a U). (4) 

For SU(2) group, a = 1,2,3. 

When the temperature is relatively high, we approximate the EOS by 

e = 3p (5) 

which is known in relativistic fluid mechanics to describe radiation. As a result, the free 
energy density can be obtained by integrating Eq. (3), 

F = ^n 4 /3 (6 ) 

where j3 is a dimensionless constant of integration. In this case, and without an external 
gluon field, the fluid system in (1) possesses a class of expanding soliton solutions which can 
be studied via variational and collective coordinate methods [10]. 

U=U iWj)> m W R ° {t T + 1)A/ " m (?) 



2 



where Rq and r are the spacial and temporal characterizations of the variational soliton and 
9{t) the usual step function in time direction. Physically, it is certainly very interesting to 
understand the origin of these two scales from a fundamental level. The approximation in 
(7) is valid provided r <C i?o- This condition enables us to define a small parameter 

For our purpose, we calculate the nonabelian current in (1) corresponding to the soliton 
solution in Eq. (7). To do so, the hedgehog ansatz is specified for the solution (7) 



U 



cos d> + ia ■ x sin t 



where x is the unit vector and (f) is given by the stereographic map 



sin (j) 



2s 



1 



COS ( 



(9) 



(10) 



1 + s 2 ' 1 + S 2 ' 

We write s as the dimensionless coordinate x/R(t). The sign in the expression of cos <fi 
signifies a topological charge which is the skyrmion number. The negative sign gives the 
skyrmion number +1 or a skyrmion and the positive sign the skyrmion number is —1 or 
an anti-skyrmion. We will take the positive sign in the following. By expressing the abelian 
current in terms of the velocity u>^ through the equation of motion, we derive the following 
expression for the nonabelian current 
/2\3 d 3 s 



d 3 xJ a » 



.pJ (1 + s 2 ) 6 
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.(11) 
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The current in (11) has a natural form of a multipole expansion due to the skyrmion orienta- 
tion in the color space. In this letter we only consider the effect of the lowest mode and the 
effects of higher polarization will be considered elsewhere. The spherically symmetric part in 
the current is contained only in the third component 

•2\3 d 3 s 



(d 3 xJ a3 ) o = -5%( 



(3 J (1 + s 2 



(12) 



where P 6 (a) = 1 - 7s 2 + 7s 4 - s 6 



3 Semiclassical gluon radiation 

Now we consider the interaction between the expanding color skyrmion and the hard partons. 
Since the transfer momentum between hard partons is in high order to that between hard 
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parton and soliton, we expect a hierarchy between the partonic coupling gyM and the effective 
coupling g e ff- Accordingly, gluon self-interaction in terms like F^ v F a ^ v can be omitted so we 
can work with a free parton picture. Then the gauge coupling in (1) becomes the coupling 
between a classical current and a free quantum field for gluon. In this approximation, the 
lowest order semiclassical amplitude is given by 

iM=g eff (l\ J d 4 xJ a »A«\0). (13) 

|0) and |1) are gluonic Fock vacuum and one-gluon state. The gluon factor in (13) is given 
by the wave function 

(l|i«(x)|0) = ^%^ (14) 

where the color and helicity parts tp, e will be summed over eventually. Putting the current 
in, we have 



iM = A(k) J dte^J {1 f S s2)6 e~ imk - s P e (s 



(15) 



where A(k) = —(2/P) 3 g e ff(p 3 £^/^/2uj. The spatial Fourier transformation can be completed 
analytically 

iM = B{k) J dte iult - R ^ k Q A {R{t)k) (16) 

where B{k) = 7T 2 A(k)/120 and Qa{x) = 5x 2 — 5x 3 + x 4 . To go further, we need to specify 
R(t) in this equation to the form given in (7). This gives 

oo 

iM = B(k)e-^ T ^ J dte^-* 4/3 Q 4 (t 4/3 ) (17) 



n 



where r\ = ujT/(kRo) 3 / 4 . With onshell condition uj = k, rj = Ak 1 / 4 where k is defined to be 
R$k. Accordingly, 

where 

oo 



(«) = J dte lAKl/4 '-' 4/3 Q 4 (t 4/3 ) (19) 



iM 

„3/4 



The radiation spectrum is given by dE = kdN '. E{k) is the total energy radiated over 
the entire time of expansion as a function of k. The number distribution is 

dM = ^\M\ 2 d 3 k (20) 

c,h 
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where the summation is over colors and helicities of the gluon. In a spherically symmetric 
setting, dAf = ndk where n is the density of states 



By straightforward calculation, 



n = 47rk 2 Y,\M^ 

c,h 



n = ai? A 2 K" 1/2 |A^ A W| 2 , 



(21) 



g = a\'^\M x (K)\\ 



(22) 
(23) 



where a = (27r 5 /225)(<7^j//3 6 ). The numerical results for A = 1/15,2/15,1/5 are given in 
Fig. 1. 
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Figure 1: Density of states and energy spectrum for A = 1/5 (Black), 2/15 (Deep Gray) and 
1/15 (Light Gray). 



4 Comparison and discussion 

Understanding the pattern of gluon radiation in relativistic heavy ion collision processes 
is important for making an accurate determination of the physical mechanisms from the 
measurement of its decay products. 

In [6], the authors extracted the asymptotic behavior of the number density in small k is 
of the \ jk form. In our case, the asymptotic of the number density in small k is ~ 1/Vk. (S ee 
Fig. 2.) The difference comes from the fact that the medium size is taken to be infinitely 
large in [6] while in our case the medium size is characterized by the soliton size Rq. So the 
IR behavior in our case is softer. 

For the case of jet quenching, the radiation energy lost is due to scattering off the hard 
quarks. A popular approach is to model the medium as a collection of colored static scattering 
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Figure 2: n/(l/ v / K) in small k for A = .2 

centers [11]. This approach can be extended to the expanding medium [4] though the gluon 
radiation by the expanding medium itself is not included. In fact, the medium induced gluon 
radiation is characterized by the frequency 



(24) 



where q is the quenching parameter, estimated to be .04 ~ .WGeV 2 / fm, and L is the in- 
medium path length of a hard parton [12]. In general loq is significantly larger than the 
characteristic momentum in our case 1/Rq. So there is a hierarchy between the medium 
induced gluon radiation spectrum and the gluon radiation spectrum by the medium. 

Our hydro dynamical approach opens up another interesting possibility to address the 
eccentricity of the elliptic flow either intrinsically by considering the nonabelian color current 
or exogenously by considering the gluon radiation patterns. This will be the topic of the 
follow-up to this work. 
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